Quantum critical behavior in the heavy Fermion single crystal Ce(Nio.935Pd .o65)2Ge 2 
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j_( ■ We have performed magnetic susceptibility, specific heat, resistivity, and inelastic neutron scatter- 

ed |' ing measurements on a single crystal of the heavy Fermion compound Ce(Nio.935Pdo.o6s)2Ge2, which 

is believed to be close to a quantum critical point (QCP) at T = 0. At lowest temperature(1.8-3.5 
K), the magnetic susceptibility behaves as x(T) - x(°) « T~ 1/6 with x(0) = 0.032 x 10~ 6 m 3 /molc 
(0.0025 emu/mole). For T < 1 K, the specific heat can be fit to the formula AC/T = 70 -T 1/2 with 
70 of order 700 mJ/mole-K 2 . The resistivity behaves as p = po + AT 3 / 2 for temperatures below 2 
K. This low temperature behavior for 7(T) and p(T) is in accord with the SCR theory of Moriya 
and Takimoto[l]. The inelastic neutron scattering spectra show a broad peak near 1.5 meV that 
?—j . appears to be independent of Q; we interpret this as Kondo scattering with Tk = 17 K. In addition, 

the scattering is enhanced near Q=(l/2, 1/2, 0) with maximum scattering at AE = 0.45 meV; we 
interpret this as scattering from antiferromagnetic fluctuations near the antiferromagnetic QCP. 

£ ■ 

In strongly correlated electron systems, a quantum critical point (QCP) separates an antiferromagnetic (AFM) or 
ferromagnetic(FM) state from a nonmagnetic Fermi liquid state at T = K. This QCP can be tuned by adjusting a 
control parameter such as doping parameter x, external pressure P, or applied magnetic field H . In the vicinity of 
the QCP, the critical fluctuations are quantum in nature and induce unique behavior. The nature of this quantum 
ground state phase transition poses one of the most significant challenges in condensed matter physics. Heavy Fermion 
(HF) compounds are very good candidates for studying the QCP. These compounds behave as Fermi Liquids (FL), 
with large values for the specific heat coefficient 7 = C/T and susceptibility x(0), and with the resistivity varying as 
Ap cx T 2 at low temperature. When such systems are tuned close to a QCP, fluctuations of the nearby magnetically 
ordered state affect the thermodynamic behavior and lead to Non-Fermi Liquid (NFL) behavior such as Ap cx T a 
^vq ■ with a < 2 and C/T cx ln(To/T) at low temperature [2]. 

The approach to a QCP has been attained in HF systems by doping in the tetragonal compounds Ce(Rui_ a; Rh 2 ;)2Si2 
(x=0.03)[3, 4] and Cei_ a: La a: Ru2Si2 (x=0.075)[5, 6] as well as in the orthorhombic compound CeCue-^Aua: (x=0.1)[7- 
9], by application of a magnetic field in YbRfi2Si2[10], and with no additional control parameter in the compound 
/3-YbAlB 4 [ll]. 

A conventional spin fluctuation theory has been established to explain the non-Fermi liquid behavior. This theory, 
implemented through a renormalization-group approach[12, 13] or through the self-consistent renormalization(SCR) 
method[l], has successfully explained various aspects of the non- Fermi liquid behavior. However, the experimental 
results for some systems do not follow this spin fluctuation theory. A locally critical phase transition has been invoked 
?H ' to explain the behavior of a few systems[14]. Hence, measurements in other systems are needed to understand the 
- - -' behavior near the QCP in HF materials. 

Below 2 K, the tetragonal compound CeNi2Ge2 exhibits behavior for the resistivity and specific heat that obeys 
the predictions of the SCR theory for a 3D AFM QCP[1]: Ap cx T 3 / 2 and C/T cx 7(0) - Ay/T [15, 16]. Inelastic 
neutron scattering shows two low energy features[17-19]. A broad peak at 4 meV that is only weakly Q-dependent 
corresponds to Kondo scattering with Tk — 46 K[20]. A peak at 0.7 meV that is highly Q-dependent and shows 
a maximum intensity at Qjv= (1/2 1/2 0) corresponds to scattering from antiferromagnetic fluctuations. Although 
this compound is clearly close to a QCP, when T decreases below 0.4 K or B increases above 2 T, the system enters 
the FL state and C/T shows saturation[15, 16, 21]. The compound can be brought closer to the QCP by alloying 
with Pd. According to the phase diagram proposed by Knebel et al[22], the QCP occurs in Ce(Ni2- 2; Pd a ;)2Ge2 when 
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x =0.065. Fukuhara et al found[23] that when x = 0.10, the ordering wavevector remains Qat= (1/2 1/2 0); for larger 
x it changes to Qjv= (1/2 1/2 1/6). 

Here we report measurements of the magnetic susceptibility, the specific heat for a series of magnetic fields, the 
resistivity, and the inelastic neutron scattering of a single crystal of Ce(Nio.935Pdo.o65)2Ge2 grown with 58-Ni to reduce 
the incoherent scattering. For this alloy, we report the inelastic neutron spectra for the first time. 

Single crystals were grown using the Czochralski method. The magnetization was measured in a commercial su- 
perconducting quantum interference device (SQUID) magnetometer. The specific heat measurements were performed 
in a commercial physical properties measurement system (PPMS). The electrical resistivity was also measured in the 
PPMS using the four wire method. The inelastic neutron experiment were performed on the MACS spectrometer [24] 
at the NIST center for neutron research (NCNR). 
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FIG. 1: (a): Magnetic susceptibil- 
ity of Ce(Ni .935Pdo.o65)2Ge2. The 
solid line is the high temperature 
Curie- Weiss fit. (b): Low tempera- 
ture inverse susceptibility. The solid 
line is x(T) - x(0) oc T~ 1/e . 



FIG. 2: (a): C/T in applied fields 
of B=0, 3 T, 6 T and 9 T. (b): Low 
temperature zero field C/T curve 
in a logarithmic temperature scale, 
(c) Low temperature zero field C/T 
data. The solid line is 70 — aT 1//2 . 



FIG. 3: Resistivity for 

Ce(Ni .935Pdo.o65)2Ge 2 . The 
solid line is the linear fit for the 
temperature range 0.4 K to 12 K 
while the dashed line represents 
T 3 / 2 dependence in the temperature 
range 0.4 K to 2 K. 



The magnetic susceptibility x(T) is shown in figure 1(a). At high temperatures the data follow Curie- Weiss (CW) 
behavior X (T) = C eff /(T - 6); in the range 150-300 K, the CW fit yields a moment 2.84 p B and 6 = -118 K. Below 
100 K, x{T) is enhanced over the CW value but no long range magnetic order is observed down to 1.8 K. In the 
traditional spin fluctuation theory [1, 2, 12, 13] for a three dimensional system near a QCP, x(T) oc T~ 3 / 2 is expected 
for an AFM while %(T) oc T~ 4 / 3 is expected for a FM. In figure 1 (b), our data for Ce(Nio.93 5 Pdo.o65)2Ge 2 follow 
the behavior of (x(T) — x(0)) oc T -1 / 6 . Exponents for the susceptibility that are smaller than 1 are also observed for 
/3-YbAlB 4 [ll] where x(T) oc T" 1 / 3 when B=0.05 T and for YbRh 2 (Sio.9 5 Geo.o5)2 where ( X (T) - X {0)) oc T-°- 6 [25]. 

In figure 2(a) we plot the specific heat C/T. The zero field data are very similar to that reported by Kuwai et al for 
Ce(Nio.936Pdo.o64)2Ge2[26]. The data are logarithmic with temperature in the range 1 K to 5 K. When an external 
field is applied, the data deviate from the InT behavior and show saturation. As for CeNi2Ge2, this is because the 
magnetic field forces the system to enter the FL state. However, in figure 2(b), the low temperature (T <1 K) data 
can be seen to deviate slightly from the InT behavior; in this temperature range the data can be fit(figure 2(c)) to the 
form 70 — aT 1 / 2 (with 70 of order 695 mJ/mole-K 2 ), which as mentioned is the expected behavior in the SCR theory. 

In figure 3(a) the resistivity p{T) data are seen to be roughly linear in temperature over a wide range 0.4 K to 
12 K, indicating NFL behavior. On the expanded scale of figure 3(b), p(T) deviates from T-linear, following a T 3 / 2 
dependence. Again, this is the expected behavior in the SCR theory. 

Hence, at temperatures higher than 2 K, this alloy exhibits typical NFL behavior, with the resistivity varying as 
p(T) oc T, the specific heat varying as C/T ocln-y?-, and the susceptibility varying as (x(T) — x(0)) oc T -1 / 6 . At 

temperatures below 1 K, however, p(T) oc T 3 / 2 and C/T oc 70 — a^/T, thereby following the expected behavior of the 
SCR theory for a three dimensional AFM QCP spin fluctuation system. 

In figure 4, we plot the INS spectra for Ce(Ni .936Pdo. 065)2662- A broad peak near 1.5 meV appears to be 
only weakly Q-dependent. As for the 4 meV peak in CeNi2Ge2, we identify this as representing Kondo scattering. 
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FIG. 4: Inelastic neutron scattering spec- 
tra for Ce(Nio.935Pdo.o65)2Ge2. The data 
were collected on MACS. In (a), the scat- 
tering at the critical wavevector Qjv = 
(1/2 1/2 0) is compared to the scattering 
at (3/4 3/4 0). The lines represent fits to 
the sum of a quasi-elastic and an inelastic 
Lorentzian peak, as described in the text. 
In (b) the difference between the data for 
Q= (1/2 1/2 0) and (3/4 3/4 0) is com- 
pared to a quasielastic Lorentzian (solid 
line) . 



Extra intensity is observed when AE < 1 meV in the vicinity ofQ = (1/2 1/2 0). This is the wavevector of the 
antiferromagnetic fluctuations seen for CeNi2Ge2 near 0.7 meV, hence we identify it as the critical wavevector Qtv 
for the QCP. 

We first fit these spectra to the sum of a quasi-elastic Lorentzian (peak width Tq) and an inelastic Lorentzian (peak 
position Ex and width Ti) (Fig. 4(a)). For Q = (1/2 1/2 0) we find £7 X =1.51 meV, Fi=0.90 meV and r =0.35 meV; 
for Q=(3/4 3/4 0), we find i?i = 1.35 meV and Fi=0.88 meV. In figure 4(b), we use a second approach to obtain the 
AFM fluctuation spectra: we subtract the spectra for Q = (3/4 3/4 0) from that measured at (1/2 1/2 0). A fit of 
the difference to a quasielastic Lorentzian peak gives the AFM fluctuation energy scale as To = 0.44 meV. 

The value of E\ is around 1.35-1.51 meV represents a Kondo temperature in the range 15.7-17.5 K. We note that 
this is smaller than the value 46 K seen in CeNi2Ge2, so that the approach to the QCP involves a reduction of Tk- 
We note in addition that Tk does not vanish at the QCP. A finite Kondo temperature is expected for a spin density 
wave system at a QCP. 

The value of the AFM fluctuation energy Tq is around 0.35-0.44 meV appears to be smaller than the fluctuation 
energy scale of pure CeNi2Ge2 where T=0.7 meV [17]. This represents "critical slowing down" of the AFM spin 
fluctutaions. However, despite the fact that this system is believed to sit at the QCP, the lifetime of the AFM 
spin fluctuations is finite, i.e. it does not diverge as expected at the critical point. Such a "saturation" of the spin 
fluctuation lifetime has been observed for other HF systems where the QCP is attained by alloying [4, 6] and may 
arise from the fact that the QCP occurs in a disordered environment. 
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